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In this work we present a partially new method to analyze 
fluctuations which are induced by causal scaling seeds. We 
show that the power spectra due to this kind of seed pertur- 
bations are determined by five analytic functions, which we 
determine numerically for a special example. We put forward 
the view that, even if recent work disfavors the models with 
cosmic strings and global 0(4) texture, causal scaling seed 
perturbations merit a more thorough and general analysis, 
which we initiate in this paper. 



I. INTRODUCTION 

At present, two ideas to explain the origin of large scale 
structure in the universe are primarily investigated: 

• Perturbations in the energy density of matter 
may have emerged from quantum fluctuations, 
which, during a phase of inflationary expansion, 
are stretched beyond the Hubble scale and 'freeze 
in' as classical fluctuations in the energy density 
of the cosmic fluid JlJ . These small initial fluctua- 
tions then evolve according to linear cosmological 
perturbation theory. 

• A phase transition in the early universe may lead 
to topological defects which seed fluctuations in the 
cosmic fluid In this case, the fluid fluctua- 
tions vanish initially and evolve according to in- 
homogeneous cosmic perturbation equations. The 
stress energy of the topological defects plays the 
role of the source term. In order for the gravita- 
tional field of the defects to be sufficiently strong 
to seed cosmic structure, we have to require e = 
AttGij 2 ~ 10 -5 , where 77 denotes the energy scale of 
the phase transition. This yields 77 ~ 10 16 GeV, a 
typical GUT scale. Topological defects which nei- 
ther over-close the universe nor die out are either 
cosmic strings or global defects. 

The anisotropies in the cosmic microwave background 
(CMB) provide an important tool to discriminate be- 
tween different models. On very large angular scales both 
classes of models lead to a Harrison Zel'dovich spectrum 
of fluctuations. For inflationary models this can easily 
be derived analytically. For defect models, the spectra 
were originally found numerically ||-J5| . Analytical argu- 
ments for this behavior are given in || . On intermediate 



scales, inflationary models predict a series of peaks due to 
acoustic oscillations in the baryon/photon fluid prior to 
recombination Q . Present observations seem to confirm 
these peaks even though the error bars are still consider- 
able [§. 

Recently, several investigations led to the conclusion 
that cosmic strings |||l0) and global O(N) defects 
do not reproduce the acoustic peaks indicated by present 
data. This led [[[2) and |l(J] to the conclusion that models 
of cosmic structure formation with scaling causal defects 
are ruled out. 

However, in a very simple parameterization of two fam- 
ilies of more general scaling causal seed models, we were 
able to fit present data very well [Q. We are thus con- 
vinced that it is too early yet to completely abandon 
scaling causal seeds as a mechanism for structure forma- 
tion. In contrary, we think it is extremely useful to study 
them in full generality, ignoring in a first step the phys- 
ical origin of the seeds. This purely phenomenological 
point of view is analogous to inflationary models where 
one sometimes manufactures the inflationary potential to 
yield the required spectrum of initial perturbations. 

To determine the power spectrum of the radiation and 
matter perturbations induced by seeds, we need to know 
the two point correlation functions of the seed energy 
momentum tensor. In this paper we present a simple 
parameterization of these functions and point out an er- 
ror frequently made in the literature. We then exemplify 
our findings with the large N limit of global scalar fields 

HQ. 

For simplicity, we work in a spatially flat Friedman 
universe. The metric is thus given by 

ds 2 = a{tf{dt 2 - rt, ,r/.r' -/./••' i , 

where t denotes conformal time. 

II. CORRELATION FUNCTIONS OF CAUSAL 
SCALING SEEDS 

We define seeds to be any non-uniformly distributed 
form of energy, which contributes only a small fraction 
to the total energy density of the universe and which 
interacts with the cosmic fluid only gravitationally. 

We parameterize the energy momentum tensor of the 
seed by 

T (seed) = M 2 Q ^ ^ (1) 
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where M denotes a typical energy of the seed and 0^„ 
is a random variable. We assume that ensemble averages 
and and spatial averaging coincide, the usual ergodicity 
hypothesis. Furthermore, we assume the random process 
M „ to be spatially homogeneous and isotropic, so that 
the two point correlation function 

(e |t „(x,t)e v (x + y,t / )) = 

i J 6 M „(x, i)0 Ap (x + y, t')d 3 x ee C^xpiY, t, t') (2) 

is a function of the difference y only. Due to the ex- 
pansion of the universe, which breaks time translation 
symmetry, we however expect C to depend on t and t' 
and not just on the difference t — t'. We consider causal 
seeds. Causality requires 



C vliXp (y,t) = 0, if |y| >t + t' . 



(3) 



The two point function in position space thus has com- 
pact support which implies that its Fourier transform is 
analytic. 

We define a seed to be scaling, if the Fourier transform, 



(4) 



contains no dimensional parameter other than t, t' and 
fcR This implies that the Ricci curvature induced by the 
seeds is a function of k and t only, multiplied by the 
dimensionless parameter e = AttGM 2 . Since we define 
Fourier transforms with the normalization 



/(k) = -J= / ./•V/.xS,x l »(,'kxS 



/(x) 



(2tt) s 



d 3 kf(k) exp(-ikx) 



(5) 
(6) 



and since M1/ (x, t) has the dimension of (length) -2 , C 
has the dimension of an inverse length. From scaling we 
therefore conclude that for purely dimensional reasons, 
we can write the correlations functions in the form 
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CV„ A p(k,M') = —F^ Xp {Vtf -k,t'/t) , 



(7) 



where F^xp is a dimensionless function of the four vari- 
ables Zi ee y/Fiki and r ee t' /t, which is analytic in zi. 

We also require the seed to decay inside the horizon, 
which implies 



1 We neglect the transition from a radiation to a matter dom- 
inated universe, which actually breaks scaling, since the scale 
t eq , the time of equal matter and radiation density enters the 
problem. In numerical examples, we have found that this 
transition in general leads to somewhat different decay laws 
for the correlation functions at large kt, but it will not alter 
our main conclusions 



C M i/Ap(k, t, t) —fkt^oo— C^ u xp(k, t, t) ^kt'^oo= . (8) 

Furthermore, since the seeds interact with the cosmic 
fluid only gravitationally, satisfies covariant energy mo- 
mentum conservation, 



(9) 



With the help of these four equations, we can, for exam- 
ple, express the temporal components, ©o^ in terms of 
the spatial ones, 0^ . The seed correlations are thus fully 
determined by the spatial correlation functions Ciji m . 
Statistical isotropy, scaling and symmetry in i, j and I, m 
as well as under the transformation i, j; k; t — > l,m;—k;t' 
require the following form for the spatial correlation func- 
tions: 

Cijlm (k, £, t ) = 

-[z i ZjZiz m F 1 (z 2 ,r) + 



Vtt' 

(zizidjm + ZiZ m 5ji + ZjZi5 lm + ZjZ m 5ti)F 2 {z 2 , r) + 
ZiZj5 lm F 3 (z 2 , r)/r + ZiZ m S ij F 3 (z' 2 , l/r)r + 
+6 ij S lm F i (z 2 , r) + (SuSjm + S im Sji)F 5 (z 2 , r)] , (10) 

where the functions F a are analytical functions of z 2 ee 
tt'k 2 , and for a ^ 3 they are invariant under the transfor- 
mation r — > 1/r, F a (z 2 ,r) = F a (z 2 . 1/r). The positivity 
of the power spectra Cijij (k, t, t) 
series of positivity conditions for the functions F a : 



0y| 2 ) leads to a 



0<F 5 (z 2 ,l), (11) 

0<F 4 (z 2 ,l) + 2F 5 (z 2 ,l) , (12) 

0< z 2 F 2 (z 2 ,l)+ F 5 (z 2 ,l) , (13) 

< z 4 Fi(z 2 , 1) + 4z 2 F 2 (z 2 , 1) + 3F 5 (z 2 , 1) , (14) 
< z^F^z 2 , 1) + 2z 2 (F 3 (z 2 , 1) + 2F 2 (z 2 , 1)) 

+F 4 (z 2 ,l) + 2F 5 (z 2 ,l) . (15) 

Since Ciji m is the Fourier transform of a real function, 



CijlmO^^t, t ) — Cijlm{ t ) 



(16) 



and thus the ansatz ( ^0[ ) implies that the functions 
F a (z 2 ,r) are real. Furthermore, decay inside the hori- 
zon (condition (|J)) yields 

lim F a {z 2 ,r)= lim F a {z 2 ,r)=0. (17) 

z 2 r — >oo z 2 /r — >oo 

In addition, analyticity implies that the functions F a do 
not diverge in the limit z — > 0, thus 

lim F a (z,r) = A a (r) 

z^O 

with 

A a (r) =A a (l/r) for all a £ 3 . 
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As an example, we have worked out the functions i*i 
to Fg in the large N limit of global scalar field seeds. 
A discussion of this simple model of scaling causal seeds 
ands its relation to the texture model of structure forma- 
tion can be found in [E]] and |l5| . In Figs. 1 and 2 we 
present the functions F^(z 2 ,r) and F2(z 2 ,r). 
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FIG. 1. The function Fs(z,r), linear (top) and \Fs(z,r)\ 
logarithmic are shown. Because of their small amplitude, the 
oscillations are virtually invisible in the linear plot. To show 
the symmetry r — > l/r, the r-axis is chosen logarithmic in 
both plots. 

The symmetry under the transition r — ► l/r is clearly 
visible. Also the conditions that F a — > if either z — > 
00 or r -t or r -» co which follows from Eq. (jlT]) is 
evidently satisfied. For fixed z the functions oscillate with 
a frequency which grows with z. Since the amplitude 
decays rapidly, these oscillations are only visible in the 
log-plots. The correlations always decay like power laws, 
never exponentially. 

The equal time correlation functions, Fi(z 2 ,l) to 
F 5 (z 2 , 1) are plotted in Fig. 3. In Fig. 4 we show A a (r). 

All the functions, except F§ which is constrained by 
Eq. (11) go through (For F\ the passage through is 
not visible since it occurs only at z ~ 30). We have also 




FIG. 2. The same as Fig. 1 for the function F2(z,r). 

verified the positivity constraints Eq. ( |l2] ) to Eq. (|lE|). 
The asymptotic behavior of the functions can be obtained 
analytically. The same is true for the functions A% to A$. 
As argued above, all the functions A4 except A3 are sym- 
metric under r — > l/r. The techniques employed to cal- 
culate the functions Fi are analogous to those explained 
in g|. 



III. SCALAR, VECTOR AND TENSOR 
DECOMPOSITION AND CMB ANISOTROPIES 

The energy momentum tensor of seeds is often split 
into scalar, vector and tensor perturbations, since the 
time evolution of each of these components is indepen- 
dent. Furthermore, due to statistical isotropy, the scalar, 
vector and tensor modes are uncorrelated. 

A suitable parameterization of this decomposition is 

k 2 1 

@ij = Sijfp ~ (kikj + 2^ Wi ^^ W j^i) Ti J ' 

(18) 

where f p , / ff are arbitrary random functions of k; w 
is a transverse vector, w • k = 0; and is a symmet- 
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FIG. 3. The functions \Fi(z, 1)| are shown. The zeros are 
visible as spikes in the log-plot. (Further below, at z ~ 30, 
also F\ passes through zero.) 



0.01 0.1 1 10 100 

r = t'/t 

FIG. 4. The functions \At(r)\ = \Fi(Q, r)| are shown. As 
discussed in the text, all of them except Az are symmetrical 
under the transition r — » 1/r. 



ric, traceless, transverse tensor, r? = Tjjfc 1 ' = 0. The 
variables /. , w and Ty represent the scalar, vector and 
tensor degrees of freedom of . The functions Fi to F 5 
determine the correlations: To work out the correlation 
functions we use 



fv 

T~ij 

where 



1 



1 



<*,,H" , 



1 



6' 
3 1 3 

3 

~2k 



Sij)Qij 



(p tl p jm - (i/2)p l3 p lm )p ma p lb e ab 

b , 

k l = k l /k 



1 ij 



(19) 
(20) 
(21) 

(22) 
(23) 



P i is the projection operator onto the space orthogonal 
to k and Sij, V { lm and ab are the projection operators 
to the scalar vector and tensor parts of Qg. 

Using these identities and our ansatz ([l0|), one easily 
verifies 



(f p (t)w*(t')) = (U(t)w*(t'))=Q 

(fp(tK j (t')) = (f p i(tH(t')) = o 

(w i (t)r jl (t'))=0 

(um(t')) 



(24) 
(25) 
(26) 
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■ ; [2F 5 (z 2 ,r)+3F 4 (z 2 ,r) 



3V«' 

+z 2 (F 3 {z 2 ,r)/r + F 3 {z 2 ,l/r)r) + 
^F 2 (z 2 ,r) + ±z i F 1 (z 2 ,r)} (27) 



1 

ftUk A 



[3F 5 + 4^ 2 F 2 + z i F 1 ] 



(28) 



-VW^F^Z 2 ^) + ^F 2 {z 2 ,r)+rF 3 (z 2 , 1/r)] (29) 
(w t (t)w*(t')) = 



- ; [F 5 + z 2 F 2 ](k 2 S l3 ~k l k 3 



k^Vtt 7 

= F 5 [Sil5j m + S im Sjl - SijSlm + k~ 2 (S l jklk m + 



^Ira^i^j fiilkjkm Si m kikj <5j/A^& m <5j"m&/A^) ~\~ 



(30) 



k kikjkikm] 



(31) 



It is interesting to note that although C'iji m is analytic, 
the correlation functions of the scalar vector and tensor 
components, in general, are not. The reason for that is 
that the projection operators S, V and T are not analytic. 
This is important. It implies, e.g., that the anisotropic 
stresses in general have a white noise and not a fc 4 spec- 
trum as erroneously concluded in p6[-^8[. The scalar 
anisotropic stress potential thus diverges on large scales, 
(I/tt| 2 ) oc l/(ifc 4 ) for kt <C 1. A result which we also have 
obtained in the large- N limit and in numerical simula- 
tions of O(N) models. The power spectrum of the scalar 
anisotropic stress potential /„• is analytic if and only if 
vector and tensor perturbations are absent, F 5 = F 2 = 0. 
In the generic situation, F$(z = 0,r = 1) = ^5(1) 7^ 0. 
We thus expect the following relation between scalar vec- 
tor and tensor perturbations of the gravitational field on 
super- horizon scales, x = kt <C 1: (The equations for 
the scalar, vector and tensor gravitational potentials in 
terms of /., w and r,j can be found in pi and pC(|.) 
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1$ - #1 



12e 2 
tk 4 
16e 2 t 



Mi) 

4s(l) 



11 ' fc 2 
|/i y | 2 }~4e 2 i 3 A 5 (l) 



(32) 

(33) 
(34) 

where \& and $ are the Bardeen potentials, <x is the vector 
contribution to the shear of the t =const. slices and 
are tensor perturbations of the metric. 

If it would be solely super horizon perturbations which 
induce the large scale CMB anisotropics, this could be 
translated into a ratio between the scalar, vector and 
tensor contributions to the Ct's on large scales, £ ~ 
50. However, since the main contribution to the CMB 
anisotropics is induced at horizon crossing, x — 1 (see 
below) the above relations cannot be translated directly 
and we can just learn that one expects, in general, con- 
tributions of the same order of magnitude from scalar, 
vector and tensor perturbations. 

Finally, we want to discuss in some detail the CMB 
anisotropics induced from scalar perturbations. In this 
case, the gravitational perturbation equations (see e.g. 
H@) imply 

$ + * = -2eU ■ (35) 

Especially, if $ has a white noise spectrum due to 'com- 
pensation' pi] , this leads to a k~ 4 spectrum for VP and 
for the combination $ — which enters in Eq. ( |36| ) . 

This finding is in contradiction with Jll|[l7j , which pre- 
dict a white noise spectrum for but it is not in con- 
flict with the Harrison Zel'dovich spectrum of CMB fluc- 
tuations which has been obtained numerically in |}] 
This can be seen by the following simple argument: Since 
topological defects decay inside the horizon, the Bardeen 
potentials on sub-horizon scales are dominated by the 
contribution from dark matter and thus roughly con- 
stant. The integrated Sachs Wolfe term then contributes 
only up to horizon scales. Therefore, using the fact that 
for defect models D g and V are much smaller than the 
Bardeen potentials on super- horizon scales (see plj), we 
obtain 

($ - X dec )jl(x - Xdec) 



(AT/T) e (k)\ sw 
+ ($' - V')(k, x)j e {x Q - x)dx , 



(36) 



where x — kt and prime stands for derivative w.r.t. x. 
The lower boundary of the integrated term roughly can- 
cels the ordinary Sachs Wolfe contribution and the upper 
boundary leads, to 



(\(AT/T)i(k)\ 2 )\ 



sir 



[3F 5 (l) + 4F 2 (l) + F 3 (l)b- 2 (x ) , 



(37) 



a Harrison-Zel'dovich spectrum. The main ingredients 
for this result are the decay of the sources inside the 
horizon as well as scaling, the rest follows for purely di- 
mensional reasons. 



IV. CONCLUSIONS 

In this paper, we outline a procedure to investigate 
causal scaling seed perturbations. We show, that the 
large number of seed correlations, which determine fully 
the induced power spectra of dark matter and CMB pho- 
tons, can be cast in only five analytic functions with cer- 
tain well defined properties. We schematically estimate 
the large scale CMB anisotropies induced. However, we 
are convinced that the relative amplitudes of large scale 
CMB anisotropies and the acoustic peaks as well as the 
dark matter power spectrum depend on details of the 
model and thus scaling causal seeds cannot be ruled out 
from studies of specific models. This finding is also con- 
firmed in |l5| ]. 

Our work is just the beginning of a program to be 
carried out which goes in different directions and to which 
we invite researchers in the field to participate. Some of 
the questions which should be explored arc the following: 

• Are there further general restrictions for the cor- 
relation functions which have not been mentioned 
here? 

• Given the functions F\ to F$ what is the exact ex- 
pression for the C/s and the dark matter power 
spectrum? What are good approximations? 

• Are there simple conditions which the functions F\ 
to F5 have to satisfy in order to lead to power spec- 
tra which are in agreement with data. 

• Are there physically plausible causal scaling seed 
models other than topological defects? 
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